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ABSTRACT 

We have selected and analysed the properties of a sample of 2905 Kg < 21.5 
galaxies in ~ 131 arcmin^ of the Great Observatories Origins Deep Survey (GOODS) 
Chandra Deep Field South (CDFS), to obtain further constraints on the evolution 
of i^s-selected galaxies with respect to the results already obtained in previous stud- 
ies. We made use of the public deep multiwavelength imaging from the optical B 
through the infrared (IR) 4.5 /im bands, in conjunction with available spectroscopic 
and C0MB017 data in the CDFS, to construct an optimised redshift catalogue for 
our galaxy sample. We computed the i^s-band LF and determined that its charac- 
teristic magnitude has a substantial brightening and a decreasing total density from 
z = to (z) — 2.5. We also analysed the colours and number density evolution of 
galaxies with different stellar masses. Within our sample, and in contrast to what is 
observed for less massive systems, the vast majority 85 — 90%) of the most mas- 
sive (M > 2.5 X 10^^ -^^o) local galaxies appear to be in place before redshift z ~ 1. 
Around 65 — 70% of the total assemble between redshifts z = 1 and z = 3 and most of 
them display extremely red colours, suggesting that plausible star formation in these 
very massive systems should mainly proceed in obscured, short-timescalc bursts. The 
remaining fraction (up to ^ 20%) could be in place at even higher redshifts z = 3 — 4, 
pushing the first epoch of formation of massive galaxies beyond the limits of current 
near-IR surveys. 
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1 INTRODUCTION 

Deep multiwavelength surveys are progressively sheding 
light on the history of galaxy evolution since very high red- 
shifts. From new surveys in the ultraviolet regime (e.g. Mar- 
tin et al. 2005) to the latest IR, sub-millimetre and radio 
campaigns (e.g. Werner et al. 2004; Dunlop 2005; Condon 
et al. 2003), the advent of numerous datasets is enhanc- 
ing our understanding of when galaxies were formed and 
how they evolved through cosmic time. UV/optical obser- 
vations can unveil the sources of stellar emission but are 
Hmited to detect unobscured systems, missing an important 
fraction of the galaxies with on-going star formation. IR 
and longer-wavelength surveys provide an unbiased way of 
studying star- forming systems, but are insensitive to galax- 
ies with little or no dust content. Near-IR observations, on 
the contrary, appear as one of the most suitable methods 
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for making a complete census of galaxy populations, as they 
are sensitive to the stellar emission from both young and old 
galaxies, and are also relatively unaffected by the presence 
of dust. The rest-frame K band is the most efficient tracer 
of the stellar mass content of the Universe. 

Recent studies of near-IR-selected galaxies allowed for a 
rapid progress in setting different constraints on the history 
of galaxy formation (Cimatti et al. 2002; Daddi et al. 2003; 
Caputi et al. 2004; Caputi et al. 2005a; among others) and, 
in particular, on the evolution of the stellar mass content of 
the Universe with redshift (Dickinson et al. 2003; Fontana et 
al. 2004; Glazebrook et al. 2004; Caputi et al. 2005a; Drory 
et al. 2005). It is now recognised that, although more than a 
half of the stellar mass has been formed at redshifts z < 1.5, 
near-IR-selected galaxies are already present before at least 
redshifts z ~ 2 — 3 and contain a significant fraction of the 
stellar mass observed at redshift z = 0. However, one of 
the main problems of existing near-IR surveys is that, when 
sufficiently deep, they are usually restricted to small areas of 
the sky, making the results prone to large error bars due to 
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cosmic variance. Also, the lack of statistics prevents refined 
studies of galaxy evolution, where trends for different near- 
IR sub-populations can be individually investigated. 

The GOODS project (Giavalisco et al. 2004) is a pub- 
lic program of multiwavelength observations which is now 
close to completion. The data sets are kindly released to 
the astronomical community in a fully reduced mode. As 
part of the GOODS project, 160arcmin^ of the CDFS 
have been observed with the Advanced Camera for Sur- 
veys (ACS) on board the Hubble Space Telescope (HST) 
in four passbands: B, V, I775 and zgso- In addition, the 
GOODS project includes deep near-IR imaging taken with 
the Infrared Spectrometer and Array Camera (ISAAC) on 
the 'Antu' Very Large Telescope (Antu-VLT). A first re- 
lease of the European Southern Observatory (ESO) Imaging 
Survey (GOODS/EIS vO.5 release) provided data in the J, 
H and Ks bands for ~ SOarcmin^ of the GOODS/CDFS 
field. In this sub-region, Roche, Dunlop & Almaini (2003) 
selected a deep sample of Extremely Red Galaxies (ERGs) 
with Ks < 22 (Vega), whose redshift distribution and lu- 
minosity evolution were studied by Caputi et al. (2004). A 
detailed comparison of the properties of ERGs with other 
-fsTs-selected galaxies was presented by Caputi et al. (2005a). 
Over the last year, however, new datasets have been made 
public within the GOODS project, increasing the coverage 
in wavelength and area of the GOODS/CDFS. The second 
GOODS/EIS (vl.O) release extends the deep J and Ks band 
imaging to ~ 131 arcmin^ of the GOODS/CDFS. Also, the 
Infrared Array Camera (IRAC) on board the Spitzer Space 
Telescope imaged this field at 3.6 /xm, 4.5 ^m, 5.8 /im and 
8.0 jj,m.. The availability of these data allows one to map the 
rest-frame T^s-band light up to redshifts z ~ 2 — 3. 

In this work, we take advantage of the further wave- 
length and area coverage of the GOODS/CDFS, in conjunc- 
tion with the information provided by the latest released 
spectroscopic and C0MB017 data, to set tighter constraints 
on the properties and evolution of Jfa-selected galaxies. The 
layout of this paper is as follows. First, in Section |21 we 
present the sample selection and give details on the multi- 
wavelength photometry. In Section |3 we explain the tech- 
nique we applied to construct an optimised redshift cata- 
logue for our galaxy sample and show the resulting red- 
shift distribution. In Section 21 we present the rest-frame 
Ks-hand luminosity function (LF) for our sample of Ks- 
selected galaxies. In Section |S] we discuss several aspects 
of the evolution of different mass systems. In Section |S| we 
explore additional hints from the study of the purely near- 
IR colours of our galaxies. Finally, in Section |7| we present 
some concluding remarks. We adopt throughout a cosmol- 
ogy with Ho = 70kms"^Mpc"\ Qm = 0.3 and Ha = 0.7. A 
Salpeter (1955) IMF over stellar masses M = (0.1-100) Mq 
is assumed, unless where otherwise stated. 



2 THE SAMPLE - MULTIWAVELENGTH 
PHOTOMETRY 

We used the public deep ISAAC Ks-hand images in the 
GOODS/EIS vl.O release to select a catalogue of Ks < 21.5 
(Vega) sources in 131 arcmin^ of the GOODS/CDFS. This 
new release is composed of 21 maps which cover more than 
twice the area of the GOODS/EIS vO.5 release. In addition. 



further improvements have been applied to the vl.O imaging, 
such as a revised data reduction technique and a uniform 
photometric calibration of the maps across the entire field. 

We performed the source extraction using the public 
code SEXTRACTOR (Berlin & Arnouts, 1996) on the Ks- 
hand images. We cut the resulting catalogue at a magnitude 
Ks = 21.5, a limit half a magnitude brighter than the cut we 
used to select our Ks sample in Caputi et al. (2005a). Sev- 
eral reasons justify our new magnitude limit. First, a few of 
the additional maps in the vl.O release have lower exposure 
times than the average exposure times of the maps compos- 
ing the vO.5 release. Second, this new limit constitutes an 
adequate balance between having access to galaxies up to 
very high redshifts and requiring for our catalogue to have 
a high level of completeness. The completeness levels of our 
vl.O sample are very similar to those estimated for the vO.5 
sample by Caputi et al. (2005a): 100%, 90% and 86% for 
Ks < 20.0, 20.0 < Ks < 21.0 and 21.0 < Ks < 21.5, re- 
spectively. Finally, our new sample has been selected hav- 
ing in mind the possibility of a complete quantitative study 
of galaxy morphology, only possible for high signal-to-noise 
objects. The results of the latter will be presented in a com- 
panion paper (McLure et al. 2005, in preparation). Our total 
vl.O Ks < 21.5 sample consists of 3184 sources. 

We measured multiwavelength photometry for our Ks- 
selected galaxies in order to model their spectral energy dis- 
tribution (SED), as we shall explain in Section |H] We ran 
SEXTRACTOR in 'double-image mode' to perform aper- 
ture photometry on the GOODS/EIS vl.O J-band images, 
centred at the position of the A's-band extracted sources. 
We measured magnitudes in an aperture size of 2.83-arcsec 
diameter for all the objects, to coincide with one of the aper- 
ture sizes used to construct the public GOODS/CDFS ACS 
catalogues. The latter have been constructed running SEX- 
TRACTOR on the zgso-band images, and making 'double- 
image mode' runs to perform photometry in the B, V and 
I775 bands. We looked for counterparts of our Ks < 21.5 
(Vega) sources in the ACS catalogues and adopted the corre- 
sponding 2.83-arcsec diameter aperture magnitudes in each 
of the four ACS passbands. A 2.83-arcsec diameter aperture 
is sufficient to collect all the fiux of the vast majority (>90%) 
of our galaxies in the different optical/near-IR passbands. 
Even in the Ks band, virtually all of our galaxies have 2.83- 
arcsec-diameter aperture magnitudes in agreement with the 
corresponding total Ks magnitudes within the error bars. 

Deep Spitzer/IRAC imaging for the GOODS/CDFS is 
now also publicly available, covering the wavelength range 
3.6 fj,m to 8.0 ^m. At the time of writing, only the first epoch 
of IRAC GOODS/CDFS data had been released. This first 
epoch includes 3.6 /xm and 5.8 /xm images for ~ 2/3 of the 
area of the GOODS/CDFS and 4.5 ^m and 8.0 /im images 
for another ~ 2/3 of this area, in such a way that only 
~ 1/3 of the field have data in the four IRAC passbands. We 
ran SEXTRACTOR on the 3.6 pm and 4.5 fim images and 
looked for counterparts of the Tfs-band sources in a matching 
radius of 1.2 arcsec. We restricted the photometry to only 
the first two IRAC channels, as the templates we used for 
the galaxy SED fitting lose accuracy or do not cover wave- 
lengths beyond 5 /im — 6 pm. We measured magnitudes in 
2.83-arcsec diameter apertures. However, the point-spread 
function (PSF) of the IRAC images is much larger than the 
characteristic PSF of the ISAAC and ACS images: the mean 
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Figure 1. The (/yrs — ^gso) versus (J — Kg) colours for the 
sources in the Ks < 21.5 sample. Cross-like symbols corre- 
spond to those objects with SEXTRACTOR stellarity parameter 
CLASS.STAR>0.8. 

full-width half maximum (FWHM) of the IRAC 3.6 ^m PSF 
is 1.66 arcsec (Fazio et al. 2004), in comparison with 0.4 arc- 
sec and 0.1 arcsec for the Ks and zsso bands in the ISAAC 
and ACS instruments, respectively. Thus, a 2.83-arcsec di- 
ameter aperture is insufficient to encompass the magnitudes 
measured in the IRAC bands with the magnitudes obtained 
in equivalent apertures on the ACS and ISAAC images. We 
computed aperture corrections for the magnitudes measured 
on the IRAC images studying the curve of growth of isolated 
stars in the field. We determined aperture corrections of 0.50 
and 0.55 magnitudes for the 2.83-arcsec-diameter aperture 
magnitudes in the 3.6 pm and the 4.5 fim bands, respectively. 

To clean our catalogues for stars and active galaxies, 
we applied criteria somewhat more strict than those im- 
plemented by Caputi et al. (2005a). We cleaned our cata- 
logues for stars using the SEXTRACTOR stellarity parame- 
ter 'CLASS-STAR', as measured on the ACS zsso images. To 
separate active galaxies, we looked for counterparts of our 
Ks < 21.5 sources in the public X-ray catalogues available 
for the CDFS (Szokoly et al. 2004). 

222/3184 sources of our initial catalogue had 
CLASS_STAR>0.8. Figure[Tlshows the location of these 222 
sources in a colour-color diagram, in comparison to the re- 
maining sources in the catalogue. We see that the major- 
ity of the objects with CLASS^TAR>0.8 occupy a clearly 
separated region in this plot, typically displaying colours 
(J — Ks) < 1.1 (Vega). However, we find 6 outliers in the 
colour-colour diagram, which we studied on a case-by-case 
basis. 4 out of these 6 outliers are X-ray classified as active 
galactic nuclei (AGN) or quasi-stellar objects (QSO). The 
remaining 2 are not detected in the X-rays, but their SEDs 
cannot be satisfactorily fitted by HYPERZ (Bolzonella, Mi- 
ralles & Pello, 2000) with any of the templates in the GIS- 
SEL98 library (Bruzual & Chariot 1993), so they are likely 



to be also stars or AGN/QSO. As a consequence, we ex- 
cluded the 222 sources with CLASS_STAR>0.8 from our 
sample. 

On the other hand, we found that 142/3184 of our Ks 
sources were detected in the X-rays. However, in principle, 
the fact that a galaxy is classified as active in the X-rays 
does not necessarily imply that its optical-to-near-IR SED is 
dominated by an active nucleus. To determine which X-ray- 
detected Ks sources had to be excluded from our catalogue, 
we modelled their SEDs using the templates in the GIS- 
SEL98 library of HYPERZ. The GISSEL98 library contains 
a wide range of synthetic SEDs based on stellar spectra but 
no power-law SED template. We determined that the SEDs 
of 57/142 X-ray-detected Ks sources could not be satisfac- 
torily fitted by any template in the GISSEL98 library, so we 
excluded these 57 objects from our catalogue, as it is likely 
that their optical/near-IR hght is at least partially contam- 
inated by an active nucleus component. 5 additional X-ray- 
detected Ks < 21.5 sources did not have either a satisfactory 
fit of their SEDs, but they already had been excluded on the 
basis of the stellarity parameter criterion. In summary, our 
final Ks < 21.5 source catalogue consisted of 2905 galaxies. 

Other no-X-ray-detected active galaxies could still be 
present in our final galaxy catalog. Another characteristic 
allowing to identify active galaxies is the power-law shape 
of their SEDs from the UV/optical through the IR regimes. 
Thus, to investigate the presence of additional contami- 
nating sources within our sample, we cross-correlated our 
galaxies with the catalogue of 24 /im-selected IR power- 
law sources in the CDFS constructed by Alonso-Herrero et 
al. (2005). Within our sample, we found only 8 no-X-ray- 
detected Ks < 21.5 sources included in the Alonso-Herrero 
et al. catalogue. We conclude, then, that the fraction of re- 
maining plausibly-contaminating sources within our sample 
is basically negligible. 



3 AN OPTIMISED REDSHIFT CATALOGUE 
FOR Ks-SELECTED GALAXIES IN THE 
GOODS/CDFS 

We used the pubhc code HYPERZ (BolzoneUa, Miralles & 
Pello, 2000) with the GISSEL98 template library (Bruzual 
& Chariot 1993) to model the SED and obtain a photometric 
redshift (zphot) for each galaxy in our sample, using photom- 
etry in up to eight broad-bands {BVIrrszssoJKs and IRAC 
3.6 j-im and 4.5 /im bands). Only models with solar metallic- 
ity have been considered. Dust-corrections have been taken 
into account convolving the template SED with the Calzetti 
et al. (2000) reddening law. 

As it was shown in Caputi et al. (2005a), the code HY- 
PERZ has a tendency to favour high- redshift estimates when 
degenerate solutions in redshift space exist. To control the 
HYPERZ output, we applied the following set of constraints 
on the obtained photometric redshifts: 

• We used the K — z relation for radio galaxies (see Fig- 
ure [IJ to determine the maximum possible redshift for each 
object. Radio galaxies delimit the behaviour of the most 
massive galaxies formed at high redshifts (e.g. Willott et al. 
2003) and normal galaxies are not usually found below this 
limit (Caputi et al. 2004; Caputi et al. 2005a). 
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Figure 2. a) Photometric versus spectroscopic redshifts for Kg < 
21.5 galaxies in the GOODS/CDFS. b) Distribution of relative 
errors. The median of this distribution is dz = (zspec — ^phot) /(1 + 
Zspec) = 0.01 and the rms is cr = 0.03. 



• The estimated redshift was constrained to the best-fit 
solution with a maximum value Zphot ~ 2 for all those galax- 
ies with a counterpart in the public shallow U-band cata- 
logues existing for the CDFS^, as higher redshift sources are 



^ Although (7-band images exist for the CDFS (Arnouts et al. 
2001), they are shallower and poorer in resolution than the ACS 
and ISAAC images used here and, thus, we decided not to in- 
corporate these data for the input catalogue of the photometric 
redshift algorithms. 
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Figure 3. The redshift distribution of Kg < 21.5 normal galaxies 
in the GOODS/CDFS. 



unlikely to be bright at such short wavelengths. Similarly, 
the estimated redshift was constrained to a maximum value 
Zphot = 4 for those sources detected in the B band. 

• We used the public code BPZ (Bem'tez 2000) to obtain 
a second, independent, set of redshift estimates for all the 
sources in our catalog. In the cases of sources with HYPERZ 
photometric redshifts Zphot > 2 not confirmed by BPZ, we 
adopted the lower estimates given by the BPZ code. 

• We used public COMBO 17 photometric redshifts (Wolf 
et al. 2004) to replace the estimates of those sources with 
magnitude R < 23.5 (Vega) and zcn < 1 (the best regime 
for COMBO 17 redshifts). 

Figure shows the resulting estimated redshifts ver- 
sus spectroscopic redshifts for all those Ks < 21.5 sources 
included in different publicly available spectroscopy samples 
in the GOODS/CDFS (e.g. Cimatti et al. 2002b; Le Fevre et 
al. 2004; Vanzella et al. 2005). We observe a very good agree- 
ment between estimated and real redshifts in most cases. 
Figure ^p) shows the corresponding distribution of relative 
errors: the median is dz — {zspec — Zphot)/{l + Zspec) = 0.01 
and the standard deviation is cr = 0.03. 

However, in spite of the good accuracy obtained for the 
photometric redshifts, the incorporation of existing spectro- 
scopic redshifts would be of benefit for the determination of 
our galaxy redshift distribution. Consequently, in order to 
optimise the quality of our final redshift catalog, we used the 
available spectroscopic redshifts for the CDFS to replace the 
photometric estimates of our galaxies whenever possible. 

The final redshift catalogue for our 2905 Kg < 
21.5 galaxies is composed of 686/2905 (23.6%) spectro- 
scopic redshifts, 629/2905 (21.7%) C0MB017 redshifts and 
1590/2905 (54.7%) HYPERZ/BPZ redshifts. 

Figure |21 shows the resulting redshift distribution of 
our sample of 2905 Kg < 21.5 normal galaxies in the 
GOODS/CDFS. 81% of these galaxies lie at redshifts z < 
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Figure 4. Observed Kg magnitudes versus redshifts z. The solid 
line corresponds to the empirical K — z relation for radio galaxies 
as obtained by Willott et al. (2003), which approximately co- 
incides with the passive evolution of a 3L* starburst formed at 
redshift z/ = 10. The dotted line is a nominal extrapolation of 
the same law. 



Figure 5. The distribution of the differences between the rest- 
frame absolute i^s-band magnitudes M^^ obtained from the clos- 
est available rest-frame band and the ones obtained extrapolating 
from the observed Kg magnitudes in all cases. The distribution 
is centered at Mk, (closest) —Mk^ (^*) = 0.04 and the rms is 
(7 = 0.24. 



1.5, while the remaining 19% is found to be at 2 > 1.5. The 
two peaks of the distribution in the redshift bins [0.5; 0.75] 
and [1.0; 1.25] are produced by known large-scale structure 
in the CDFS (e.g. Le Fevre et al. 2004). The resuhing Hub- 
ble diagram Kg versus z is shown in Figure |3] 



4 THE EVOLUTION OF THE i^s-BAND 
LUMINOSITIES 

Caputi et al. (2005a) computed and analysed the rest-frame 
Kg-hanA luminosity function (LF) from redshifts (z) = 1 
to {z) = 2.5. In this work, we take advantage of the larger 
area covered by the i^s-band GOODS/EIS vl.O release to 
set tighter constraints on the galaxy rest-frame Kg-hanA LF, 
and to model its evolution up to (z) = 2.5. The incorpora- 
tion of IRAC 3.6 /im and 4.5 data enables us to compute 
this LF with no need of magnitude extrapolations up to 
redshifts z ~ 0.6 and z ^ 1, respectively, and with smaller 
extrapolations at higher redshifts. Anyway, we note that the 
k-corrections in the Kg-haxiA are less dramatic than those 
in other shorter wavelength bands (Poggianti et al. 1997). 

We obtained the rest-frame absolute magnitude Mk^ of 
each galaxy directly from the HYPERZ output. For those 
objects with spectroscopic, C0MB017 or BPZ redshifts, 
we adopted the best SED fitted by HYPERZ constrained 
to the predetermined redshift. For the remaining objects, 
the best-SED fit and redshift estimation were simultane- 
ously obtained. HYPERZ computes the rest-frame absolute 
magnitude in the selected filter from the observed passband 
which is the closest to that filter in the rest-frame. In our 
case, the rest-frame Mk^ magnitudes are obtained either 



from the Kg, 3.6 /^m or 4.5 /im apparent magnitudes, de- 
pending on the redshift of the source and the availability 
of IRAC data. With this 'closest band' method, magnitude 
extrapolations are minimised. 

It is interesting to compare the rest-frame absolute Kg- 
band magnitudes Mk^ obtained with this procedure, and 
those obtained when no IRAC data is taken into account 
for the SED fitting, i.e when the rest-frame Kg-hanA mag- 
nitude extrapolations are directly made from the observed 
Kg-haxiA in all cases (as, for instance, in Caputi et al. 2005a, 
where IRAC data were not available) . Figure |S] shows the 
distribution of the difference between the absolute magni- 
tudes Mks obtained with the two different methods. The 
statistics has been made over the magnitude range rele- 
vant for the computation of our LF (see below), i.e. over 
all the sources with —27 < Mk^ < —22, as obtained with 
the 'closest-band' method used here. We see that, although 
the resulting distribution has a long tail, the standard devi- 
ation a is small: the center of the distribution is located at 
Mk, (closest) -Mk, (Ks) = 0.04 and the rms is cr = 0.24. 
Thus, we conclude that, in most cases, the rest-frame abso- 
lute magnitudes Mk, obtained extrapolating the observed 
Kg-ha.nA magnitudes are consistent with those derived from 
the longer-wavelength data used here. 

We computed the rest-frame Kg-haxiA LF using two dif- 
ferent techniques: the l/Vmax method and a maximum like- 
lihood (ML) analysis (MarshaU et al. 1983). In the 1/V,nax 
method, the LF is directly computed from the data, with 
no parameter dependence or model assumption. The ML 
analysis, on the other hand, is a parametric technique which 
requires an a priori functional form for the LF and its evolu- 
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Table 1. The best-fit parameters for the ML analysis of the rest-frame A'a-band LF. The values of a and M*{z = 0) have been taken 
from Kochanek et al. (2001). a is considered to be constant with redshift. The remaining parameters have been taken as free for the 
modelling, ipoiz = 0) is the LF normalisation parameter, for which we obtain a best- fit value in excellent agreement with Kochanek et 
al. (2001). 



a 


M*(z = 0) 




ku 


ZM (/.o(2 = 0) (Mpc-^) 


-1.09 (constant) 


-24.16 ± 0.05 


1.67 ±0.08 1.78 ±0.06 


0.63 ± 0.10 


1.88 ±0.10 (3.9 ± 0.3) X 10-3 



tion with redshift. We parametrised the rest-frame ifs-band 
LF with a Schechter function (Schechter 1976): 



L 



L 



LHz) 



0.4{M-Af (z)) 



exp 



L*{z) 



L*{z)) ■■ 
(1) 



and we assumed that the density and luminosity evolved 
with redshift as 



t>o(z) — 4>o{z = 0) X exp 



M*(z) = M*(2 = 



V ZM ) 



(2) 



(3) 



The parameters a and M*{z — 0) have been fixed to the 
local value: a = -1.09 and M* {z = 0) = -24.16 ± 0.05 
(Kochanek et al. 2001). The obtained best-fit values for the 
remaining (free) parameters are given in Table 

In Figure|S| we present the resulting rest-frame A's-band 
LF at redshifts (2) = 1.0, 1.5, 2.0 and 2.5. The empty sym- 
bols indicate the LF obtained with the 1/Vmax method. The 
values shown are strictly above the luminosity complete- 
ness limits of each redshift bin. The number of sources in 
the luminosity-completeness range are 822, 335, 145 and 83 
for the redshift bins [0.75; 1.25], [1.25; 1.75], [1.75; 2.25] and 
[2.00; 3.00], respectively. To correct for the small incomplete- 
ness of our sample to Ks = 21.5, we weighted each source 
with a multiplicative factor depending on its apparent mag- 
nitude: 1.00, 1.11 and 1.16 for sources with Ks < 20.0, 
20.0 < Ks < 21.0 and 21.0 < Ks < 21.5, respectively (see 
Section |5J. These correction factors have a minor impact on 
the resulting LF (cf. Caputi et al. 2005a). The error bars 
correspond to the maximum of Poissonian errors and the 
errors due to cosmic variance. We set a fixed value of 15% 
of the number counts for the cosmic variance error bar at 
z > 1 (cf. Somerville et al. 2004). The solid lines show the 
evolution of the LF with redshift, as obtained through the 
ML analysis. In the four panels of Figure |H] we have added 
for comparison the local K-band LF measured from Two 
Micron All Sky Survey (2MASS) data by Kochanek et al. 
(2001; dotted lines). The LFs obtained by the two different 
methods are in good agreement in the different redshift bins. 

From the upper- left panel in Figure |S| we can confirm 
once more that the bright end (Mx^ < — 25) of the rest- 
frame ifs-band LF significantly increases from redshifts z = 
to 2 = 1. This trend has already been observed in previous 
works (Drory et al. 2003; Feulner et al. 2003; Pozzetti et al. 
2003; Caputi et al. 2004; Caputi et al. 2005a). From our ML 
analysis we find that, although the total density of objects 



(j)o decreases by a factor ~ 1.5 between redshifts z — and 
z = 1, the characteristic magnitude M* brightens in ~ 0.7 
magnitudes between those redshifts. The overall result is a 
rise in the density of objects with Mk, < — 25. 

Recently, Dahlen et al. (2005) found that the charac- 
teristic magnitude M* of the rest-frame J-band LF suffered 
a dimming from redshifts z ~ 0.4 to z ~ 0.9. This result 
was apparently in contradiction with the brightening found 
by Feulner et al. (2003) for the same LF in a similar red- 
shift range. However, as Dahlen et al. (2005) explain, the 
one-to-one comparison of the different parameters charac- 
terising the Schechter function is misleading, as M* , a and 
(po are coupled, and different behaviours in the evolution of 
M* can still produce consistent LFs. Interestingly, on the 
other hand, Dahlen et al. (2005) infer that, given the dim- 
ming they find for M* in the rest-frame J-band LF, the rest- 
frame J^s-band LF should present a similar behaviour. They 
suggested that the brightening observed by different authors 
in the rest-frame A's-band LF might be due to inaccuracies 
in the k-corrections used to extrapolate from the observed 
to the rest-frame Ks-hand at z ~ 1, rather than being a 
real effect. The incorporation of IRAC 3.6 /^m and 4.5 /xm 
data in the present work allows us to determine that this is 
not the case. Once more, one should conclude that a one- 
to-one comparison of the best-fit parameters obtained for 
the Schechter function (and even more if they are inferred) 
is misleading. The comparison of the densities of objects in 
each luminosity and redshift bin found by different surveys 
is the correct way of concluding on the resulting LFs. 

As it was shown by Caputi et al. (2004; 2005a), after 
the rise of the bright end of the rest-frame Ks-hand LF from 
redshifts z = to (z) = 1, the very bright end {Mk^ ^ ~ 
26) of this LF does not show any sign of decline up to at 
least redshift {z) — 2.5. Our ML analysis suggests that the 
density of the Mk^ J; ~ 26 objects might increase up to 
(z) = 1.5 and then stay constant up to (2) — 2.5 within the 
error bars (cf. Figureinj. Again, this can be explained by the 
combined effect of the number and luminosity evolution we 
obtain through the ML analysis: the total density of objects 
00 decreases by a factor ~ 6 between 2 = to (2) = 2.5, 
but the characteristic magnitude M*{z) brightens by more 
than one magnitude between these redshifts. 

The combined density and luminosity evolution we de- 
termined through the ML analysis also implies a particularly 
strong decrement in the number of intermediate and low lu- 
minosity {Mks > —25) objects with increasing redshift. We 
remark, however, that we do not directly see these inter- 
mediate and low luminosity objects at high redshifts, given 
the limits of our survey. We only extrapolate this behaviour 
from the LF obtained through the ML analysis, which as- 
sumes the validity of a Schechter function for the description 
of the LF at all redshifts. Thus, two alternative conclusions 
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Figure 6. The rest-frame Kg-band luminosity function at redshifts (z) = 1.0, 1.5, 2.0 and 2.5 (circles, squares, up triangles and 
diamonds, respectively). The empty symbols correspond to the LF computed with the l/Vmax method. Only values above the luminosity 
completeness limits of each redshift bin are shown. The solid lines show the resulting LF obtained through the ML analysis. For a 
comparison, we show in each panel the local K-band LF measured from 2MASS data by Kochanek et al. (2001; dotted lines). 



can be extracted from this result: 1) a significant fraction of 
low and intermediate luminosity objects are created between 
(z) = 2.5 and « = 0, or 2) at least as many low luminos- 
ity objects arc present at high as at low redshifts, and our 
assumptions for the computation of the LF with the ML 
analysis are not valid to describe the faint-end of the LF 
at high redshifts. The direct access to low luminosity ob- 
jects at high redshifts through deeper surveys is necessary 
to conclusively prove the behaviour of the faint-end of the 
LF. 



In hierarchical models, the creation of new objects is 
the consequence of the production of galaxy mergers and 
subsequent triggering of star formation. Very recently, Con- 
sclicc (2005) found that the fraction of mergers in faint and 
low mass galaxies peaked at lower redshifts than for higher 
luminosity and more massive galaxies. This result appears 
to be consistent with the presence of luminous galaxies at 
high redshifts z ^ 2.5, as they would be produced from 
mergers at those or oven higher redshifts. The decline in the 
number of intermediate luminosity objects at high redshifts 
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4.0x10'm, <M< 1.5x lo'" M, 1.5 x lo'" M, < M < 3.0 x 1o"'m, 3.0xU)"'m < M < 6.0 x lo'" M 




Figure 7. The observed (/tts — Ks) colours versus redshift for galaxies within different ranges of stellar mass. In all the panels, circles 
and triangles indicate precise and lower limits to the (7775 — Ks) colours, respectively. The horizontal dashed line represents the typical 
extremely-red-galaxy colour cut, {7775 — Ks) = 4, while the vertical dotted line delimits the region of completeness for each stellar mass 
bin. 



that we extrapolate from the ML analysis would also be 
consistent with the fact that a significant fraction of them 
form through mergers at later times. However, the differen- 
tial redshift peak observed for mergers of different mass and 
luminosity still cannot decide on the formation epoch of all 
the lowest luminosity objects, and the two hypotheses dis- 
cussed above could still be possible. At redshifts z <1, ma- 
jor mergers become progressively less important (e.g. Patton 
et al. 2002), so it is expected that the evolution of the Ks- 
band LF is mainly governed by the dimming of the already 
present old stellar populations. The significant decrease of 
the bright end {Mr, < -25) of the LF from redshifts {z) = 1 
to 2 = is probably indicating that the most massive galax- 
ies present at 2 ~ 1 must have finished the construction of 
the bulk of their stars by this redshift, and mainly fade out 
all the way through redshift z = Q (cf. Section [5.21 . 

5 THE EVOLUTION OF THE STELLAR MASS 

5.1 Computation of stellar masses 

We used the rest-frame Ks-haxiA luminosity of each galaxy 
to compute its stellar mass. We obtained the rest-frame ab- 
solute magnitude Mk^ (and corresponding luminosity) from 
the IfYPERZ best-fit SED on an individual basis, as ex- 



plained in Section 0| We used the public code GALAXEV 
(Bruzual & Chariot 2003) to construct grids of mass-to- 
light [M/Lks) ratios depending on galaxy age, assuming a 
Salpeter IMF over stellar masses M = (0.1 - 100) M©. The 
advantage of using luminosities in the rest-fame Ks-hanA is 
that other parameters, such as the dust-correction or the 
galaxy star-formation history, have a minor importance in 
the determination of the M/ Lk^ ratios. As in Caputi et 
al. (2005a), we modelled the evolution of these ratios using 
two galaxy templates, with no dust, corresponding to dif- 
ferent exponentially declining star-formation histories with 
characteristic times r ~ 0.1 Gyr and r ~ 5Gyr. We deter- 
mined the M/ Lks of each galaxy using the model most simi- 
lar to the HYPERZ best-fit template in each case, and inter- 
polating between the values which corresponded to the clos- 
est ages in the grids. Our estimated stellar masses are typi- 
cally accurate within a factor < 2 (where the larger uncer- 
tainties correspond to plausible degeneracies in parameter 
space, especially in those cases in which the galaxy redshift 
is simultaneously determined with the best-SED fit). On the 
contrary, the dependence on star-formation history and dust 
corrections would be much stronger if optical rather than 
near-IR luminosities were used to derive the stellar masses. 

As we discussed in Section |3] the rest-frame Ks-hanA 
luminosity derived as here and the luminosity that would 
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be derived from the extrapolation of the corresponding 
observed Ks-hanA magnitude are consistent within a rms 
a = 0.24. This implies that the derived stellar masses will be 
consistent within IQ"''-*^ (=1=0.24) ^ j g within a factor ~ 1.25. 
Thus, we conclude that the use of the observed B through 
T^s-band data, with no incorporation of IRAC data, would 
have introduced an additional uncertainty of ~ 25% on the 
derived stellar masses. This is a relatively minor contribution 
to the error budget, which is mainly dominated by plausi- 
ble degeneracies in parameter space and especially by the 
adoption of a fixed IMF. 



5.2 The colours and stellar masses of galaxies at 
different redshifts 

Figure 13 shows the observed (7775 — Ks) colours versus red- 
shift for galaxies with stellar masses within different mass 
ranges. In all the panels, circles represent those objects with 
precise (7775 — Kg) colours, while triangles indicate lower 
limits to the colours of objects with magnitude 7775 > 26 
(Vega; equivalent to < 2a detections in the 7775 band), com- 
puted as (26 — Ks). The horizontal dashed line shows the 
typical colour cut defining extremely red galaxies (ERGs; 
[hib — Ks) > 4, Vega), while the vertical dotted lines delimit 
the region of completeness for the lowest masses in each mass 
bin. The mass completeness limits have been estimated us- 
ing the median of the k-corrections and the maximal M/ Lk^ 
ratios at each redshift (corresponding to a galaxy formed in 
a single burst at redshift zj = 00). At redshifts z > 1, typi- 
cal ERG colours indicate the presence of old stellar popula- 
tions and/or young stellar populations heavily enshrouded 
by dust. On the contrary, {1775 — Ks) < 4 colours ai z > 1 im- 
ply the presence of unextincted on-going star-formation. At 
z < 1, the observed (7775 — Ks) — 4 colour cut is merely for- 
mal, as only systems with very large dust extinctions would 
show colours above that limit. 

Inspection of Figure|7|shows that the (hrs—Ks) colours 
of galaxies at a given redshift depend on the stellar mass. 
The fraction of red-to-blue galaxies systematically increases 
with assembled stellar mass, even considering each mass bin 
only up to the completeness redshift limit. In particular, the 
vast majority of galaxies with M > 2.5 x 10^^ Mq between 
redshifts z — 1 and z = 4 are ERGs. Some of these galax- 
ies might be purely passive systems (cf. e.g. Daddi et al. 
2005), while in other cases the extremely red colours would 
only be explained by the superposition of old stellar compo- 
nents and dust (see Caputi et al. 2004; Papovich et al. 2005; 
and cf. Section |SJl. The extremely red colours characteris- 
ing these very massive galaxies reveal important constraints 
on their star-formation histories. In contrast to what is ob- 
served in other less massive systems, our results indicate that 
unobscured time-extended star-formation at z > 1 is rare in 
galaxies with assembled stellar mass M > 2.5 x 10^^ Mq. 
The main build-up of these systems must be produced in 
short bursts and/or with considerable amounts of dust. 

Figure |5] shows an amplified version of the lower-right 
panel of Figure[7| For reference, we added in this figure the 
modelled evolution of the (7775 — Ks) colours for galaxies 
formed in an instantaneous burst at redshift Zf = 5 and 
Zf = 3 (solid and dotted-dashed lines, respectively), both 
with no dust and passive evolution thereafter. From Figure|S| 




Figure 8. The observed (7775 — Kg) colours versus redshift for 
galaxies with estimated stellar mass M > 2.5 X 10^1 Mq. The 
symbols are the same as in Fig. |7| The solid and dotted-dashed 
lines show the modelled evolution of the (7775 — Kg) colours for 
galaxies formed in an instantaneous burst at redshift zj = 5 and 
Zf = 3, respectively, both with no dust and passive evolution 
thereafter. 



we can see that, at 2; > 1, galaxies within the ERG regime 
with assembled stellar masses M > 2.5 x lO^'^ Mq display a 
wide range of (7775 — Ks) colours. This suggests that they 
must undergo a variety of star-formation histories, either 
different formation redshifts or at least different formation 
redshifts for the youngest stellar population present in them. 

At redshifts z < 1, on the contrary, the galaxies with 
stellar mass M > 2.5 x 10^^ Mq present in our sample show 
quite homogeneous (7775 — Ks) colours. This is again in con- 
trast to what is observed in less massive systems. The ho- 
mogeneous colours displayed by these most massive galax;ies 
suggest that their optical/near-IR light must be dominated 
by old stellar populations. However, this does not necessarily 
mean that all of them are following a strictly passive evo- 
lution, but that any possible on-going star formation would 
only produce minor changes to their properties (e.g. stel- 
lar mass). In fact, a substantial fraction ( > 30%) of these 
galaxies are detected in the mid-IR, i.e. experience on-going 
star- formation (Caputi et al. 2005b, 2005c). The typical IR- 
derived star-formation rates of these mid-IR detected galax- 
ies at z < 1 are 10 - SOMeyr"^ (Caputi et al. 2005b). If 
we consider that the lifetime of a starburst episode is at 
most ~ 10® yr, we infer that these galaxies should be form- 
ing 1 — 5 X 10^ Mq per starburst episode at those redshifts, 
i.e. an additional minor amount to their already existing 
stellar mass. Some of these very massive systems at 2: < 1 
could be candidates for cluster E/SO galaxies, as early-type 
cluster galaxies show very little dispersion in their colours 
from z = (e.g. Bower et al. 1992) to 2; ~ 1 (e.g. Ellis 
et al. 1997; Stanford et al. 1998; Blakeslee et al. 2003). In- 
deed, the existence of large-scale structure is suggested by 
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Figure 9. The comoving number density of galaxies with stellar 
mass M > 1 X lO^^M© (circles) and M > 2.5 X IO^^Mq (triangles) 
(Salpeter 1955 IMF, restricted to M = (0.1-100) Mq). The filled 
triangle at z = 3.5 denotes the re-computed number density of 
M > 2.5 X 10^"*^ Mq galaxies after correcting for the effects of large- 
scale structure in the corresponding redshift bin. The star-like 
symbols indicate the density values for the same mass thresholds 
at 2 = 0, obtained integrating the local stellar mass function (Cole 
et al. 2001; Bell et al. 2003). The dashed lines show the evolution 
of the number density of haloes with total mass M > 2 X 10^^ Mq 
and M > 1 X IO^'^Mq (upper and lower curves, respectively), as 
obtained from A-CDM models of structure formation (with a code 
kindly provided by Will Percival). The dotted line corresponds to 
the predicted number density of galaxies with mass M > 1 X 
10"'"'^ Mq obtained with SPH simulations (Nagamine et al. 2005). 

the colour-redshift diagram, as a considerable number of our 
M > 2.5 X 10^^ Mq galaxies appear to lie at discrete redshift 
values. On the contrary, only a morphology study would be 
able to conclude on the plausible early or late-type classifi- 
cation of these systems (McLure et al. 2005, in preparation). 

5.3 The evolution of the galaxy number and mass 
densities 

As we did in Caputi et al. (2005a), we analyse the evolu- 
tion of the number density of different-mass galaxies over 
cosmic time. In the present work, however, due to the anal- 
ysis of a considerably larger area, we are able to select 
a statistically significant sample of very massive galaxies 
(M > 2.5 X IO^^Mq), which enables us to explore the very- 
high mass end of the stellar-mass evolution. 

Figure|5|shows the comoving number density of galaxies 
with assembled stellar mass M > 1 x 10^^ Mq (circles) and 
M > 2.5 X IO^^Mq (empty triangles). The asterisks denote 
the corresponding local values, as obtained from the inte- 
gration of the local stellar mass function (Cole et al. 2001; 
Bell et al. 2003). As in the computation of the galaxy LF 
(cf Sectional, the error bars correspond to the maximum 
of Poissonian errors and the errors due to cosmic variance. 



Figure 10. The comoving number density of galaxies with stellar 
mass M > 1 X IO^^Mq (Chabrier 2003 IMF). The symbols are 
the same as in Figure 1^ 

All the results presented in this section incorporate Vmax 
corrections, but deliberately not sample incompleteness cor- 
rections (which anyway are very small). We decided to do 
so to explicitly show that the contrast betweeen some of 
the observational and theoretical results discussed below is 
independent of such corrections. 

The evolution of the densities of galaxies with stellar 
mass M > 1 X lO^^Afg shows that a minimum of 15% of the 
local value is in place before redshift z — 3 (this value rises to 
~ 20% considering the number densities cut at stellar masses 
M > 1.3 X 10^^ Mq, i.e. the regime of strict mass complete- 
ness to z — 3). This results is consistent with Caputi et 
al. (2005a). However, in contrast to this previous work, we 
observe here that the number density of Af > 1 x 10^^ Mq 
systems does increase from redshifts z « 3 to z 1.5, al- 
though this increment is modest (~ 20% of the local value). 
More than 50% of the local Af > 1 x 10^^ Mq galaxies appear 
to have formed at redshifts z < 1.5. 

The dashed lines in Figure |U] trace the evolution of 
dark- matter haloes with total mass M > 2 x 10^^ Mq and 
M > 1 X IO^^Mq (upper and lower curves, respectively), 
as obtained from the extended Press-Schechter formalism 
(Bond et al. 1991), with a code kindly provided to us by 
Will Percival. The mass thresholds for the halo distributions 
have been deliberately selected to coincide with the corre- 
sponding number densities of galaxies at 2 = 0, assuming 
a single halo occupation number. The compared observed 
and predicted density evolutions confirm the trend found in 
Caputi et al. (2005a): at redshifts 1 < z < 3, the density 
of haloes with total mass M > 2 x IO^'^Mq is systemati- 
cally larger than the density of galaxies with stellar mass 
M > 1 X 10" Mq. Caputi et al. (2005a) also showed that 
the number density of galaxies and haloes coincide again at 
higher redshifts z ~ 3.5, although this cannot be re-assessed 
from the present sample as the Ks — 21.5 cut implies com- 
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pleteness only for stellar masses M > 2.5 X 1O"M0 up to 

On the other hand, we investigate the very high mass 
end (M > 2.5 x 10^^ M0) of galaxy number density evo- 
lution and also compare with the corresponding theoret- 
ical predictions. Inspection of Figure |5] shows that up to 
^ 20% of the local massive galaxies with assembled stellar 
mass M > 2.5 X 1O"M0 is in place before a = 4, while 
the construction of the majority of these very massive sys- 
tems occurs between redshifts 2 = 3 and z = 1.5. Within 
our sample, we find that ~ 90% of the local galaxies with 
mass M > 2.5 x lO^^M© is already in place by redshift z ~ 1. 
This is contrast to what is observed for less massive systems, 
whose evolution is only halfway complete by that redshift. 

The comparison with theoretical models shows that, if 
we assume the same conversion between baryonic and dark- 
matter halo mass as inferred at 2 = 0, there appear to be 
too few dark matter haloes of adequate mass in place at high 
redshift to host the observed number of high mass galax- 
ies. We observe that the deficit is particularly important at 
z > 3. Similar incompatibilites between the presence of very 
massive systems and the predictions of theoretical models 
at high redshifts have already been pointed out in previous 
studies (e.g. Fontana et al. 2003; Genzel et al. 2003; Pozzetti 
et al. 2003; Cimatti et al. 2004; Fontana et al. 2004; Mc- 
Carthy et al. 2004; Saracco et al. 2005). Our results show 
that this incompatibility might certainly exist for the most 
massive galaxies. It could be argued, however, that our ex- 
cess of very massive galaxies might be biased by cosmic vari- 
ance effects (more than what is accounted for within the er- 
ror bars). In effect, as we noted in Section [5.21 the existence 
of large-scale structure within our sample can be deduced 
from the colour-redshift diagrams. Within our sample, we 
find 16 galaxies with M > 2.5 x IO^Mq at 3 < 2 < 4, yield- 
ing a number density pc = (3.87 ±0.97) x lO"'^ Mpc"^. 9/16 
of these galaxies have photometric redshifts concentrated in 
the narrow interval 3.55 < z < 3.65, and there is also a spec- 
troscopically confirmed QSO at redshift z = 3.61 within our 
surveyed field (not included in the final galaxy sample anal- 
ysed in this work) , strongly suggesting the presence of large- 
scale structure for very massive galaxies at these redshifts. 
To obtain a lower limit on the number density of galaxies 
with stellar mass M > 2.5 x IO^Mq at (2) = 3.5, we re- 
computed this density excluding the 9 galaxies at 3.55 < 
z < 3.65. We obtained = (1.69 ±0.64) x 10-^Mpc-^ i.e. 
logio(pc) = -4.77+0 21 (filled triangle in Figure l^, while 
the density of the corresponding haloes at that redshift is 



only logio(pc 



-5.33. Thus, it seems that even taking 



into account the effects of large-scale structure, we cannot 
completely reconcile the theory predictions with the obser- 
vations at high redshifts. 

Recently, Nagamine et al. (2005) showed that the prop- 
erties of high redshift red and massive galaxies could be well 
reproduced through the use of hydrodynamical simulations 
with star-formation rate histories that peak at higher red- 
shifts than those used in semianalytic models. It is then in- 
teresting to compare also their predictions with our observed 
evolution of number densities of massive galaxies. The dot- 
ted line in Figure |5] shows the predicted number density of 
galaxies with stellar mass M > 1 x 10^^ Mq, as obtained 
with the Nagamine et al. smoothed particle hydrodynamics 
(SPH) simulations. We observe that these simulations yield 
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Figure 11. The total stellar mass densities obtained from the 
present study (large circles and up-triangle, with error bars, 
for precise values and lower limit, respectively), compared to 
values previously obtained by other authors: Dickinson et al. 
2003 (down-triangles); Fontana et al. 2003 (asterisks); Rudnick 
et al. 2003 (diamonds); Glazebrook et al. 2004 (right-triangles); 
Fontana et al. 2004 (squares and double squares for the observed 
and extrapolated values in the K20 survey, respectively); Ca- 
puti et al. 2005a (small circles and double circle for total and 
M > lO^'^ Mq mass densities, respectively); Drory et al. 2005 
(crosses). The circle with a star at Zpf^^t = show the local value 
of the total stellar mass density obtained from the integration of 
the local galaxy stellar mass function (Cole et al. 2001, Bell et al. 
2003). The solid line is the evolution of the stellar mass density as 
obtained from the integration of the star-formation rates derived 
from Sloan Digital Sky Survey (SDSS) data (Heavens et al. 2004). 
Dickinson et al. values correspond to the average of the values ob- 
tained with the 1-component and 2-component models with solar 
metallicity (see Table 3 in Dickinson et al. 2003). Rudnick et al. 
values are computed only on galaxies with rest-frame V-band lu- 
minosity Lv > 1.4 X IO^^Lq and, thus, miss a significant fraction 
of the mass. Glazebrook et al. values only correspond to galaxies 
with estimated mass log(M) > 10.46 and are at least partially 
incomplete above redshift z ~ 1.2. All the values in this figure 
assume a Salpeter IMF over stellar masses M = (0.1 — 100) Mq 
and a cosmology with Ho = 70kms~^Mpc~^, Cm = 0.3 and 
= 0.7. 



densities below the observed values, and the discrepancies 
are larger at low than at high redshifts. 

However, in the Nagamine et al. models, galaxies masses 
are directly obtained from the simulations and are IMF in- 
dependent. Thus, for a proper comparison with these model 
results, we need to test whether the adoption of a differ- 
ent IMF to compute our stellar masses could improve the 
agreement between model and data. Figure 1101 shows the 
evolution of the comoving number density of galaxies with 
stellar mass M > 1 X 10" M0 considering a Chabrier (2003) 
IMF. Both our observed densities and the local value have 
been transformed in accordance to this IMF. The change 
from a Salpeter (1955) to a Chabrier (2003) IMF produces 
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an overall re-scaling of all our stellar masses by ~ 0.3 dex, 
in such a way that a cut M > 1 x 10^^ for masses in the 
Chabrier IMF corresponds to a cut M > 2 x lO^M© in the 
Salpeter IMF. 

We repeat the same exercise as in Figure |^ and com- 
pare the observed values with the densities of dark matter 
haloes with a mass threshold deliberately selected to coin- 
cide with the density of galaxies at z = 0. As in Figure |31 
we observe that, at very high redshifts z > 3, the predicted 
number of haloes appears to be insufBcient to encompass the 
observed densities of M > 1 x 10^^ M© (Chabrier) galax- 
ies. On the contrary, the SPH predicted number densities 
of M > 1 X 10" M© galaxies by Nagamine et al. (2005) ap- 
pear to be in quite good agreement with our observed values 
from redshifts z ~ 1 through 2 ~ 4, although they produce 
some excess of massive local galaxies. We note that a similar 
behaviour would be deduced from the densities of dark mat- 
ter haloes if we did not constrain the comparison to those 
haloes with densities coinciding with the local densities of 
massive galaxies. This shows that what we find here is a di- 
rect discrepancy with the predictions of hierarchical models, 
and does not depend on the specific recipes adopted either 
by semianalytic models or SPH simulations. 

In summary, hierarchical models do not seem to be able 
to reproduce the observed densities of the most massive 
galaxies from the local Universe to very high redshifts. Al- 
ternatively, it could also be argued that the use of a single 
IMF is not suitable to convert from luminosities to stellar 
masses in the whole range of redshifts 2 = to z = 4 (see 
e.g. Kroupa & Weidner 2005, for a discussion). Unluckily, 
the validity of a single IMF is an issue extremely hard to 
test at high redshifts, and it is not possible to overcome 
the IMF uncertainty when doing the comparisons between 
observations and theory. 

Finally, we computed the evolution of the total stel- 
lar mass density with redshift, as derived from our sample 
of Ka < 21.5 galaxies. Figure [TTI shows our obtained stel- 
lar mass densities (large circles and up-triangle, with er- 
ror bars), compared to other values previously published in 
the literature. From the present sample we find that the to- 
tal stellar mass densities at (z) — 1.25 and (z) = 1.75 are 
Pc = (2.34±0.35) X lO**M0Mpc-^ and pc = (1.33±0.20) x 
lO**M0Mpc-^ respectively, i.e. (48 ± 7)% and (27 ± 4)% 
of the local value. These values very well agree with the 
previous determinations made by Caputi et al. (2005a) and 
other works. At higher redshifts, the completeness limits of 
our sample only allows us to obtain lower limits on the total 
stellar mass densities. The up-triangle in Figure 1111 repre- 
sents a lower limit to the total stellar mass density at red- 
shift (z) = 2.5. As our galaxy sample is only complete for 
masses M > 1.3 x 10^^ Mq up to z ~ 3, we compare this 
lower limit to the stellar mass density contained in galaxies 
with stellar mass Af > 1 x 10^^ Mq obtained by Caputi et 
al. (2005a) (double circle in Figure ITT^ . and observe a very 
good agreement. The total stellar mass densities at redshifts 
z > 2 have been better constrained by previous deeper sur- 
veys. The estimated stellar mass densities at z ~ 2.5 and 
z ~ 3.5 are (20-24)% and (13-18)% of the local values, re- 
spectively (Caputi et al. 2005a; Drory et al. 2005). 
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Figure 12. The observed versus predicted (Ks — [3.6 /.im]) and 
{Kg — [4.5/jm]) colours (upper and lower panels, respectively) 
for the Ks < 21.5 (Vega) galaxies in the restricted area of the 
GOODS/CDFS studied by Caputi et al. (2005a). The dashed 
lines indicate the region of consistency, taking into account the 
error bars on the observed colours. 



6 ADDITIONAL HINTS FROM PURE 
NEAR-IR COLOURS 

Caputi et al. (2005a) predicted the pure near-IR Ka/IKAC 
colours of Ti's-selected galaxies based on the best-fit SED of 
each object. These predictions were made for the galaxies 
selected in a sub-region of 50.4 arcmin^ of the area analysed 
in this work. In the light of the currently available IRAC 
data, we can assess the validity of these predictions. Fig- 
ure 1121 shows the observed versus predicted {Ks — [3.6fj,m]) 
and (Ks — [4.5 fim]) colours for the Ks-hand galaxies in the 
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Figure 13. The observed {Ks - [3.6 ^tm]) and [Ks - [4.5/im]) 
colours (upper and lower panels, respectively) versus redshift for 
the Ks < 21.5 (Vega) galaxies in the GOODS/CDFS. 



restricted GOODS/CDFS area (upper and lower panels, re- 
spectively). For clarity, we only show those sources with 
small error bars in their observed IRAC magnitudes, i.e. 
e(3.6 fim) < 0.3 and e(4.5 ^m) < 0.3. As in all cases the Ks 
magnitudes have e{Ks) < 0.2, the measured colours are typi- 
cally precise within ~ 0.4 mag (dashed lines in both panels of 
Figure IT^ . From thisfi gure we determine that 89% and 84% 
of the predicted {Ks — [3.6 /^m]) and {Ks — [4.5 /^m]) colours, 
respectively, are in agreement with the observed colours, 
within the error bars. We can confirm, then, the majority of 
the near-IR colours predicted in Caputi et al. (2005a). 

Figure 1131 shows the observed {Ks — [3.6 /im]) and 
{Ks — [4.5 /im]) colours versus redshift (upper and lower 



panels, respectively) for the Ks < 21.5 (Vega) galaxies in 
the larger sample studied in this work. 3.6 and 4.5 fim 
data are available for alternative ~ 2/3 of the sources in our 
sample. Once more, for clarity, we only show those sources 
with magnitude errors e(3.6 fim) < 0.3 and e(4.5 ^im) < 0.3. 
The variation of these observed colours with redshift con- 
firms the trends predicted by Caputi et al. (2005a) within 
the error bars. Inspection of Figure [T^ reveals. however, the 
existence of a population of low-redshift {z < 0.5) galaxies 
with colours {Ks — [3.6^im]) > 1.0 and {Ks — [4.5jj,m]) > 1.5, 
which are quite redder than those predicted for low-redshift 
galaxies in the previous work. The colour excess of these 
galaxies cannot be reconciled with the predictions by simply 
taking into account the photometry error bars. The analysis 
of these sources shows that 31 out of the 32 objects sat- 
isfying z < 0.5 and {Ks — [3.6 /im]) > 1.0 and all of the 
17 objects satisfying z < 0.5 and {Ks — [4.5 /im]) > 1.5 
are very low mass systems, with assembled stellar masses 
M < 2.5 X 10^ Mq- The majority of them are best-fit by 
young-galaxy templates. The excess of near-IR colour could 
be explained by the presence of an active nucleus, whose 
IRAC magnitudes could not be properly predicted using 
the GISSEL98 stellar templates. On the other hand, an al- 
ternative and perhaps more likely explanation is that these 
very red Ks/IKAC colours are produced by the presence of 
prominent emission lines in the IRAC bands. The presence 
of an aromatic feature at 3.3 ^m is characteristic of star- 
forming galaxies. These aromatic features are not incorpo- 
rated in the GISSEL98 templates, explaining why such very 
red colours could not have been predicted from the SED 
models. If this is the case, the colour excess would be the 
signature of intense star-formation in these young objects. 
Only near-IR spectrocopy would be able to confirm the na- 
ture of these sources. 

Figure ITU shows the observed {Ks — [4.5 fim]) colours 
versus redshift in different stellar mass bins, analogously to 
figure |7| We observe that the {Ks — [4.5 /xm]) colours in- 
crease up to redshift 2 ~ 2, and the colour-redshift relation 
becomes tighter with increasing stellar mass. In contrast, 
the colour dispersion becomes quite large above that red- 
shift, where the observed near-IR bands start to map the 
rest-frame optical bands. This colour dispersion is in part 
due to larger photometry uncertainties, but it is probably 
also an indication of a wide variety of star-formation histo- 
ries and dust extinction in these galaxies. 

Figure \W\ shows an amplified version of the high-mass 
panel in figure 1141 We added for comparison the modelled 
evolution of the observed {Ks — [4.5 /^m]) colours for galaxies 
formed in a single burst at redshifts Zf = 3 and 2/ = 5, both 
with no dust and passive evolution thereafter. Some inter- 
esting conclusions can be extracted from the analysis of this 
figure. Firstly, it is quite evident that the observed near-IR 
colours of many massive galaxies at high redshifts {z > 1.5) 
cannot be explained without considerable amounts of dust 
(cf. also Caputi et al. 2004). Independently of the formation 
redshift, the colours of a passively-evolving single burst are 
not sufficient to reproduce the red observer-frame near-IR 
colours of massive galaxies. We conclude, then, that purely 
passive galaxies should be a minority at high redshifts. Sec- 
ondly, the colours observed at lower redshifts {z < 1.5) in- 
dicate that the near- IR light of M > 2.5 x 10^^ Mq galaxies 
must be dominated by evolved stellar populations. A sim- 
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Figure 14. The observed {Ks — [4.5 fj,m]) colours versus redshift for galaxies within different ranges of stellar mass. In each panel, the 
vertical dotted line delimits the region of completeness for each stellar mass bin. 



liar conclusion has been obtained from the analysis of the 
(^775 — Ks) colours of these systems. 

We note that the tight colour-redshift relation observed 
for high-mass galaxies at z < 1.5 — 2.0 does not imply a 
common formation redshift or star-formation history. For in- 
stance, as we mentioned in Section [5.21 some of these galax- 
ies are bright mid-IR sources, while others are not, suggest- 
ing a wide variety of on-going star-formation rates (Caputi 
et al. 2005b,c). However, and in contrast to what is pro- 
gressively observed in less massive systems, all the plausible 
new star-formation in very massive galaxies at low redshifts 
could have a minor impact on their overall properties and, 
in particular, on their SED colours. 



7 SUMMARY AND CONCLUDING REMARKS 

With the GOODS project close to completion, in this work 
we have been able to set tighter constraints on the evolution 
of near-IR-selected galaxies in the GOODS/CDFS field. We 
took advantage of the latest available very deep and homo- 
geneous multiwavelength photometry, in conjunction with 
the information provided by complementary spectroscopic 
programs and the C0MB017 campaign, to construct an op- 
timised redshift catalogue for the Ks < 21.5 galaxies in the 
GOODS-South field. Based on this catalogue, we performed 



a detailed analysis of galaxy luminosity and mass evolution 
up to very high redshifts z ^ 3 — 4. 

The incorporation of the Spitzer/IRAC data at 3.6 and 
4.5 fj,m allowed us to compute rest-frame Ks-hand luminosi- 
ties and stellar masses with considerably smaller extrapola- 
tions than those applied in Caputi et al. (2005a). In spite 
of this important difference, we note that most of the main 
conclusions obtained by Caputi et al. (2005a) have been con- 
firmed by this work. The addition of IRAC data did not pro- 
duce any overall change in the derived galaxy luminosities 
and stellar masses. Other authors also have concluded that 
the use of IRAC data produced a minor impact on the de- 
rived galaxy stellar masses, when other shorter optical/near- 
IR multiwavelength data are also available (e.g. Shapley et 
al. 2005). 

From the modelling of the evolution of the Ks-hand 
luminosities, we determined that the characteristic magni- 
tude M* of the rest-frame Ks-hand LF significantly bright- 
ens from redshifts z = to (z;) = 2.5. This is simulta- 
neously produced with a decrease in the total density of 
objects, so we conclude that the Ks-hand LF has both a 
substantial luminosity and density evolution with redshift. 
The characteristics of this evolution reveal some important 
aspects of the sequence of galaxy formation. For instance, 
that most of the brightest galaxies are assembled at high 
redshifts (2 > 1), while the creation of a substantial num- 
ber of intermediate-to-low luminosity objects might continue 
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Figure 15. The evolution of the observed (^Kg — [4.5 /im]) colours 
with redshift for galaxies with estimated stellar mass M > 2.5 X 
10^^ Mq. The solid and dotted-dashed lines show the modelled 
evolution of the [Ks — [4.5 /^m]) colours for galaxies formed in an 
instantaneous burst at redshift Zf = b and = 3, respectively, 
both with no dust and passive evolution thereafter. 



down to low redshifts (although the latter will need to be 
confirmed by future deeper surveys, which have direct access 
to low- luminosity objects at high redshifts). 

We have also studied the very high stellar-mass end 
(M > 2.5 X 10" Mq) of galaxy evolution. We found that 
most of these very massive systems in our sample have been 
formed in the ~ 3.5 Gyr-period of time between redshifts 
2 ~ 3 and z ~ 1, although a significant fraction (up to 
~ 20%) appear to have been in place before redshift 2 = 4. 
Very recently, from the cross-correlation of our Ks < 21.5 
galaxy sample with the Spitzer/24 /^m galaxies in the same 
field, Caputi et al. (2005b) found that the most-intense star- 
forming galaxies at redshifts z ~ 2.0 — 3.0 are mainly massive 
(M > 1.5 X lO^^M©) galaxies. This is also the main epoch 
of flourish of submillimetre galaxies (Chapman et al. 2005), 
which presumably are progenitors of the most massive local 
galaxies (e.g. Stevens et al. 2003) . Thus, this period appears 
to be critical to understand the formation and evolution of 
the bulk of the most massive systems present in the local 
Universe. 

Theoretical models of galaxy formation in their current 
status are not completely able to account for all of these 
results. Several different clues have been proposed in recent 
works to try to improve the agreement of the models with 
the data at high redshifts, including e.g. the effects of QSO 
and supernova feedback (Granato et al. 2004; Cirasuolo et 
al. 2005), the adoption of an alternative IMF (Baugh et al. 
2005) or the use of hydrodynamical models (Nagamine et 
al. 2005). However, the discrepancies are not produced by 
the specific recipes used by theoretical models to account 
for star formation, but seem to lie at the more fundamental 
level of hierarchical model predictions. Hierarchical models 



do not seem to be able to reproduce the whole of observa- 
tional properties from the local Universe to high redshifts. 
A likely alternative to this problem is the invalidity of a uni- 
versal IMF when deriving stellar masses from the observed 
luminosities. A definitive comprehension of the main epoch 
of very massive galaxy formation, as well as the determina- 
tion of the initial distribution of stars in these systems, still 
constitutes a major challenge from both the observational 
and theoretical points of view. 
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